Novel phase transition in collective motion with appearance of abnormal agents 
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We introduce a novel type of abnormal agents that proceed in the opposite direction of that defined for the 
normal agents. A new order parameter, y, is introduced to describe the characteristic of the system. Many inter- 
esting phenomenons emerge as the number of abnormal agents number shifts, e.g., the system may transform 
into a new phase (from y~ltoy~-l) suddenly with abnormal agents getting denser, or stay disordered 
forever. A variety of other properties like size effect, agent speed and sight radii that have impacts on the new 
collective dynamics are also studied in detail. We suggest that our model or its modified versions can be applied 
to explain a variety of phenomenons with multiple kinds of particles interacting with each other and shape the 
system dynamics. 
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I. INTRODUCTION 

Collective motion is a fascinating phenomenon that 
emerges in the nature from microscopic scale like molec- 
ular motors im , bacterialH^and individual cellslH H], to 
macroscopic level such as insect swarms||31 and bird flocksfl ■ 

In 1995, Vicsek et al. have proposed an elegant and 
powerful model to describe such behaviors ^. Due to its 
elegancy and efficiency, Vicsek Model (VM) has received 
a wide attention, and its mutations have been extensively 
discussed. In Ref. lllOll . impacts of alignment rules or inelastic 
collisions were studied. Ref . lUlll studied the influence 
of simultaneously presented volume exclusion and self- 
propulsion. The collective motion of polar units moving in 
two dimensions with nematic collisions was presented in 
Ref. III2I1 . In Ref . IllsUlql adhesion between the particles was 
introduced to avoid isolation of clusters under open boundary 
conditions. The optimal view angle for the formation of 
collective motion was studied in Ref lflTl [18 [. Adaptive 
speed to accelerate flocking was studied in Ref. 111911 . Ref. ll20ll 
provided a comprehensive review of the development of 
research in collective motion, from the basic ideas to exciting 
recent discoveries. 



one type of self-propelled agents, i.e., all the agents share the 
same properties. Nevertheless, it is interesting to question 
that, if we place agents with different behaviors together, 
how will their interactions shape the system dynamics. 
Introducing different types of agents in a modified VM 
Model may provide us a potential tool kit for a wider range 
of research topics. For example,the formation of mutual 
cooperative phenomenon, the spread of infectious disease 
between different species, a unique form of colony shaped by 
the interaction between two types of coexisting bacteria and 
the diffusion of atoms in a certain kind of material. Therefore 
we propose that further study of the collective dynamics with 
multiple types of self-propelled particles may both deeper our 
understanding in theoretical nonequilibrium statical physics 
and provide useful methods for practical modeling of the 
nature. 

This paper is arranged as follows. In SecHH we de- 
scribe the rules of out model in detail, and illustrate the 
states the system may be in. In SecHIll we present how 
the phase transition occur with parameters running. In 
SeclIVI we inspect further into the system dynamics with a 
number of numerical simulations. At last, a summary is given. 



In most of these modified versions, there are usually only 



II. MODEL WITH ABNORMAL AGENTS 
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A. Model description 

In our model, all the agents are assumed to move in a 
square-shaped cell of linear size L with periodic boundary 



conditions. Initially, each agent is randomly distributed in 
the square, with its velocity direction randomly distributed in 
[tt, -7t). The position of the ith agent is updated as 



Xi(t + 1) = Xi(t) + Vi{t) 



(1) 



at each time step. For a normal agents, its direction is defined 
as in VM: 



0i{t + 1) ^ {0i(t))r + A0i, 



(2) 



where (OiU)},- is the direction of the average velocity of the 
agents within the horizon radii r of the ith agent, including the 
ith agent itself. It is obtained from: 



(6»,(f)>,- = arctan[(v,sin6»,(f)>,/<v,cos0,(f)>.]- 



(3) 



In EqlJl A9i is a random number evenly distributed in [-/;, /;), 
denoting the thermal noise of the agents. Unless otherwise 
noted, T] is set as 0. In other words, we mainly focus on ana- 
lyzing noise-free systems. 

Besides the normal agents, we also introduce a new type of 
agents named the abnormal agents. The only difference be- 
tween the normal agents introduced in the VM and the abnor- 
mal agents we introduce is that, the direction of the velocity 
of an abnormal agent at time f + 1 is opposite to the direction 
of the average of the velocity of agents around it at time t: 



0i{t + 1) = -{0i{t)}r + A0i. 



(4) 



This new property of the abnormal agents is simple but 
nontrivial. The dynamic of a single abnormal agent is no 
more complicated than a normal one. However, this property 
makes the abnormal agents behave completely contrary to 
the normal ones in the same environment. We choose this 
property not only because we want to study the system with 
two type of particles behaving in a totally contrary way, but 
also because it can be a representative example for further 
studies involving abnormal agents with different properties. 

In each simulation, we place «„;, abnormal agents and «„,„ 
normal agents in the system. In the latter section we will 
demonstrate how phase transition occur with n„i (and n,„„, 
sometimes) shifting. 

We define the order parameter (p as the normalized average 
velocity: 



1 



^nm{ab) 



Nnm(ab)^'0 



ienni(ab) 



(5) 



where the subscript nm and ab stand for the normal agents and 
the abnormal agents, respectively. 

A new parameter, y, is introduced in our work to depict which 
phase the system is in. It is defined as: 



nm{ab) lenm{ab) 



ienm(ab) 



From the expression one will see that y always satisfies -1 < 
y<l. 




(a)Phase 1: n„„, = 800, n„i, = 200 




(b)Phase 2: n,„„ = 200, n^i, = 800 

FIG. 1: Phase 1 and 2 collective motion states in two consecutive 
time steps 




(a)Collective motion with n„„, = 570, «„(, = 930 




(b)Collective motion with n„„, = 530, »„(, = 970 
FIG. 2: Collective motion states near phase transition point 

B. Different phases of the system state 

With nab and n,„„ shifting, we can observe the transition be- 
tween two phases of collective motion. One of them is similar 
to the collective motion in VM, while the other is a novel type 
that has never been reported. The representative figures are 
presented in FigH] 



The first phase (referred to as Phase 7) of collective 
motion, similar to the one appears in the VM, will happen in 
our model when iiab < n„m or nab ■^ n,„„. In this phase, all 
the normal agents form a steady flow toward one direction, 
as illustrated in Fig |l(a)| in two consecutive time steps. 
Obviously y„,„ -- 1 in this phase. 

The second phase (referred to as Phase 2) happens when 
nah > n„,„ to a certain extent. This is a novel state that 
does not exist in literature. In phase 2 the abnormal agents 
outnumbered the normal ones and dominate the collective 
dynamics of the system. Fig |l(b)| illustrates this state in two 
consecutive time steps. 



In the first step, all the abnormal agents head toward the 
same direction, and the normal agents towards the op- 
posite. 

In the following step, as abnormal agents will adjust 
their direction opposite the average direction of the pre- 
vious step, which is actually the the average direction of 
the abnormal ones as they outnumber the normal ones, 
all of them turn around and head opposite to their pre- 
vious direction. 




-Abnormal agent 
- Normal agent 



n =900, n =100, L=25 



FIG. 3: Disordered state, with n,„„ = 100, Uat = 900 and L = 25 



2. So do the normal agents as they tend to follow the aver- 
age direction of the agents around each agent in the first 
step. 



III. PHASE TRANSITIONS 



A. Phase transition with respect to ;in(, 



3. In the third step, the system goes back to the state in 
step one, and this circulation goes on forever 

In this phase, as both types of agents are always turning 
around, we have y,,,,, -- -1. Therefore, y„„, can be used as a 
parameter to describe the phase transition. 

To further illustrate the system dynamics near the phase 
transition points, we present Fig[2l Fig |2(a)| illustrate a 
phase 1 state near the transition point. One can see that the 
abnormal agents in the area free from the flow of normal 
agents have formed states similar to phase 2, but locally in 
small groups. The normal agents form a cluster dense enough 
to sail through the sea of abnormal agents, which would bring 
the normal agents into phase 2 if the normal agents is not 
dense enough. This happens just after we increase iiah by 50 
and decrease «,„„ by 50. The system transformed into phase 2 
is depicted in Fig |2(b)| . 

Besides forming collective motion, the system can also 
stay disordered forever, if the densities of the abnormal agents 
Pah and the normal agents p,„„ are low enough, as illustrated 
in Fig [3] In this situation, the abnormal agents dominate 
the system dynamic, and form phase 2-like state locally. 
However, the information of the direction of velocity cannot 
spread to the globe, and the whole picture is still of no order 



As indicated in FigH] and |2] the phase transition appears 
after «„/, outnumbers «„„,. In the previous section we have 
proposed y,,,,, as a parameter to describe the phase the system 
is in. Fig m illustrates the relation between y„„, and iiat, with 
«„,„ fixed as 500. This just reflect the fact that the ratio of Uab 
and n„,„ decides which type of agents paly the dominant role 
of the collective dynamic, so the collective dynamic shows 
different patterns with respect to changing agents number 
ratio. 



B. Phase transitions with respect to size efTect 

Size effect is an important topic when discussing self- 
propelled particles. We are aware whether the system 
dynamic is directly related to L or not, and will demonstrate 
that this relation do not exist in the following. 

In Figl3] one may notice that not only the ratio, but also 
the absolute values of p„/, and p,„„ play a critical role in the 
formation of collective motion. For phase 1, we have already 
understood its similarity with the VM. So in this subsection 
we only focus on the discussion of the situation in phase 2. 
In Fig|5], we fix «„/, - 950 and «„,„ = 50, and let L run from 
3 to 23. One can find that with L increasing and density of 
agents decreasing, the system is less likely to form collective 
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FIG. 4: First order phase transition of y with growing iiab 




FIG. 6: y with respect topa/,, with different values of L. 
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FIG. 5: Times reach collective motion with respect to L. At each 
value of L 500 ?? runs have been performed 
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motion. This reflects the fact that when pahlpmn is large but 
pab too small, the phase 2 state cannot form globally. Figl3] 
gives one typical screenshot of this state. 



FIG. 7: y with respect to jiat-, with different values of v. 



we describe can emerge at any scale. 



IV. FURTHER INSPECTION INTO THE SYSTEM 
DYNAMICS 

The influence of L 

Fig|6]illustrates the phase transition with respect to pat, with 
different values of L. In the simulations, p„,„ is fixed as 6. 
One can find that with L running larger to infinity, the curve 
become stable. This proves the phase of the system to be 
determined by the ratio of p„„, and pab- (Of course, according 
to Fig|5] the absolute value of p„„, and p„i are also critical.) 
Further more, the system dynamics is not directly related to L, 
which leads to a nontrivial conclusion that the phenomenon 



The influence of velocity 

Besides L, we are also interested to see how the value of 

V can influence the system dynamics. Fig]?] illustrates how 
the critical points moves to the right side with v increasing. 
What's more, the range in which the system has a chance 
to stay in either phase 1 and phase 2 gets larger, e.g., when 

V - 0.01, there is only one point with -0.7 < y < 0.7, but 
there are three when v = 0.03. Each of these points represents 
average of 50 runs, and each run produces a y with |y| ~ 1 . 
This implies that with v vanishing, the phase transition goes 
more sharply. 

The Influence of sight radii 

FigS depicts how the critical point change with respect to 
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when abnormal agents outnumber normal agents. When 
r - inf, the phase transition happens exactly at p^b - prim. 
This is understandable, since all the agents share the same 
(0,(0) when r = inf, and (0,(0) is completely decided by the 
side with more agents. 

The influence of thermal noise 

In the previous analysis, the thermal noise 77 is not taken into 
consideration. At the end of this paper, we present Fig|9]to 
show the impact on (p„„, from rj. One can conclude from Fig|9] 
that when «,„„ and itab bring the system near the critical point, 
e.g., nab - 500 or riab - 1000, the system become more likely 
to be disordered with 77 increasing. On the contrary, when 
iiab = or Tiab = 1500, the system stays far from the critical 
points and shows in a higher capability to remain stable. 



FIG. 8: y with respect to nab, with different values of r. 
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V. SUMMARY 

We have introduced a novel type of agents in the self- 
propelled particle system. We have described the rules of 
the model in detail, and introduced the concept of abnor- 
mal agents and order parameter y. With y indicating the 
characteristic of the system, we have discovered some new 
types of phase transitions with respect to the values of n„,„ 
and iiab- The abnormal agents play the role of blocking the 
normal agents and impeding the formation of the previously 
known collective motion. However they can also promote 
the formation of another type of collective motion with 
their number growing. What's more, we have studied the 
properties of the system dynamics in detail, and presented a 
variety of numerical simulation results. We suggest that the 
model we introduce may imply a new research sub-direction: 
the interaction of agents with different behaviors in one 
system. 



FIG. 9: <p„,n with respect to t], with different values of riab- 
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emerges earlier when «„/:, increases. The curve also goes down 
more sharply with a larger r. With r running to the infinity, 
the phase transition became of first order, and happen exactly 



We would like to give our appreciations to Wen-Yao Zhang 
and Yuki Nagato for their kindly help. This work is funded 
by the National Natural Science Foundation of China (Grant 
Nos. 9102402610975126, 10635040 and 11005001). 



[1] Campas O., Kafri Y., Zeldovich K. B., Casademunt I. and 
Joanny J.-F, Phys. Rev. Lett., 97,038101(2006). 

[2] Segel L., and SIAM, I. Appl. Math, 32,653(1977). 

[3] Allison C. and Hughes C., Sci. Prog., 75,403(1991). 

[4] Shapiro I. A. and Dworkin M., Bacteria as Multicellular Organ- 
isms (Oxford University Press 1997) . 

[5] Nedelec F, Surrey T, Maggs A. and Leibler S., Nature (Lon- 
don), 389,305(1997) . 

[6] Gregor T., Fujimoto K., Masaki N. and Sawai S., Science, 



328,1021(2010). 
[7] Wilson E. O., The Insect Societies (Belknap Press of Harvard 

University, Cambridge, Mass. 1971). 
[8] Parrish I. K. and Hamner W. M., Animal Groups in Three Di- 
mensions (Cambridge University Press, Cambridge) 1997. 
[9] T. Vicsek, Phys. Rev. Lett 75, 1226(1995). 
[10] D. Grossman, I. S. Aranson, and E. B. lacob, New lournal of 

Physics 10, 023036(2008). 
[11] F Peruani, A. Deutsch, and M. Bar, Physical Reveiw E 74, 



030904(R)(2006). 
[12] F. Ginelli, F. Peruani, M. Bdv, and H. Chate, Physical Review 

Letters 104, 184502(2010). 
[13] G. Gregoire, H. Cliate, and Y. Tu, Physica D 181, 157(2003). 
[14] H. Chate, F. Ginelli, G. Gregoire, F. Peruani, and F. Raynaud, 

Eur. Phys. J. B 64, 451(2008). 
[15] P Szabo, M. Nagy, and T. Vicsek, Physical Review E 79, 

021908(2009). 
[16] L. Q. Peng, Y. Zhao and B. M. Tian. et al. Physical Review E 



79,026113(2009). 
[17] B. M. Tian, H. X. Yang and W. Li. et al Phys. Rev. E 

79(2009)052102. 
[18] Y J. Li, S. Wang and Z. L. Han. et al. EPL, 93 (201 1) 68003 
[19] J. Zhang, Y. Zhao and B. M. Tian. et al. Physica A 388 (2009) 

1237C1242. 

[20] T. Vicsek and A. Zafiris, larXiv: 1010.50171 /1]. 



